Polarization-dependent ponderomotive gradient force in a standing wave 
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The ponderomotive force is derived for a relativistic charged particle entering an electromagnetic 
standing wave with a general three-dimensional field distribution and a nonrelativistic intensity, 
using a perturbation expansion method. It is shown that the well-known ponderomotive gradient 
force expression does not hold for this situation. The modified expression is still of simple gra- 
dient form, but contains additional polarization-dependent terms. These terms arise because the 
relativistic translational velocity induces a quiver motion in the direction of the magnetic force, 
which is the direction of large field gradients. Oscillation of the Lorentz factor effectively doubles 
this magnetic contribution. The derived ponderomotive force generalizes the polarization-dependent 
electron motion in a standing wave obtained earlier [A.E. Kaplan and A.L. Pokrovsky, Phys. Rev. 
Lett. 95, 053601 (2005)]. Comparison with simulations in the case of a realistic, non-idealized, 
three-dimensional field configuration confirms the general validity of the analytical results. 
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I. INTRODUCTION 

The ponderomotive force is a time-averaged force expe- 
rienced by a charged particle in an oscillating electromag- 
netic (EM) field that is spatially inhomogeneous. In the 
standard perturbative approach [H 0], it is shown that 
a charged particle in an oscillating EM field attains an 
oscillatory quiver momentum superimposed on a slowly 
varying guiding center momentum p. The latter is sub- 
ject to the classical ponderomotive force Fp-. 



dt ~ 



(1) 



where m is the mass of the particle, e its charge, 
eo the permittivity, uj the frequency of the EM field 
and I{x) the position-dependent field intensity. The 
classical ponderomotive force is of gradient form, and 
always directed toward regions of low field intensity. 
The ponderomotive force is observed and exploited 
in a wide range of contexts. In laser-plasma physics, 
the ponderomotive force drives the formation of laser 
wakefields that are used for next generation electron 
accelerators [sl, 1^ . Ion beams are produced by intense 
laser irradiation of thin foils, in which the ponderomotive 
force plays an essential role [1, Schemes have been 
proposed for ponderomotive laser-vacuum acceleration 
of electrons ,7, 8]. In Paul traps, ions are confined 
by a ponderomotive potential [9|. In electron beam 
diagnostics, the length of electron bunches is measured 
by sequentially scattering different sections of the bunch 
using the ponderomotive force of a laser pulse [HI • 

The field gradients that can be obtained in a single 




FIG. 1: Charged particle with initial velocity vo, which is 
deflected by the ponderomotive force of a standing EM wave 
oriented with its nodal planes parallel to the [x, j/)-plane. The 
dashed arrow indicates the polarization direction for the case 
considered in Section [III Al 



laser pulse are set by the laser pulse duration longi- 
tudinally and the focal spot transversely. For many 
applications of the ponderomotive force this means 
that, in order to obtain a sufficiently strong force, 
field intensities are required that are large enough 
to cause a relativistic quiver motion (which happens 
if the normalized amplitude of the vector potential, 
a EE eA/mc = e^/2I / {tQc) / {mau) > 1, i.e. / > 2 • 10^® 
W/cm^ for a wavelength of 800 nm). Relativistic field 
intensities necessitate more complicated descriptions of 
the average EM force [l2| - [l3 |. or at least restricts the 
domain of validity of Eq. ([1]) [l^ [l^ . An intermediate 
situation occurs when an already relativistic particle 
enters an EM field with nonrelativistic intensity. A 
relativistic derivation 1 171 shows that this introduces an 
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additional factor \/l + / (mc)^ + in the denomi- 
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nator of Eq. ([T]), resulting in an accurate description for 
practical situations [Tsj . 

An alternative to the application of a relativistic laser 
pulse is the use of a standing wave. In this configuration 
the nodes and antinodes are spaced on the scale of 
the wavelength, resulting in large field gradients. For 
example, a standing wave produced by two counter- 
propagating EM waves with wavelength A = 800 nm 
and a very modest, nonrelativistic peak field intensity 
of 10^^ W/cm^, already causes ponderomotive forces of 
the order oi Fp/e ^ 1 GV/m. For this reason, a number 
of applications of the ponderomotive force have been 
proposed that take advantage of the large field gradients 
in a standing EM wave. Hebeisen et al. [l^ suggested 
a table-top standing wave version of the bunch length 
measurement setup mentioned previously. Following an 
earlier idea [l^], Balcou proposed a novel X-ray free 
electron laser based on the wiggling of electrons in the 
ponderomotive potential of a standing wave [21]. Faure 
et al. used a standing wave formed by colliding laser 
pulses to pre-accelerate electrons ponderomotively in 
a laser- Wakefield setup [l^, demonstrating that the 
production of monoenergetic electron beams can be 
made stable and reproducible in that way [23l - |25| . 
Baum and Zewail proposed to create attosecond electron 
pulse trains by bunching of an electron beam due to a 
co-moving ponderomotive beat potential between laser 
pulses of different frequencies [2g . 



ponderomotive force in a standing wave is polarization- 
dependent. Consider Fig. [1] showing a particle with 
charge q and initial velocity Do parallel to the x-axis, 
incident on a standing wave with electric field E 
and magnetic field B. The wave is oriented with its 
nodal planes parallel to the (x, 2/)-plane, so that the 
spatial variation of the field is much more rapid in the 
z-direction than in the transverse direction. This is 
the typical system considered in this paper. When the 
particle enters the EM field, it will start to quiver in 
the polarization direction in response to the oscillating 
electric force qE. This 'electric quiver', combined with 
the Lorentz force equation, leads to the well-known 
average force, Eq. ([Ij, independent of the polarization 
direction. However, the incident particle will also quiver 
in response to the magnetic force sa q{vQ x B). This 
'magnetic quiver' can be comparable in magnitude to 
the electric quiver for relativistic particles. Because 
the magnetic quiver is in the z-direction, the particle 
samples a large field gradient, leading to an additional 
contribution to the ponderomotive force that is com- 
parable to the electric one. And since the magnetic 
force and hence the amplitude of the magnetic quiver 
depend on the angle between and B, the magnetic 
contribution is dependent on the polarization direction. 



In view of all these important technological applica- 
tions, a thorough understanding of the ponderomotive 
force in a standing wave is essential. The scattering of 
charged particles by a standing EM wave was first de- 
scribed in a quantum-mechanical context by Kapitza and 
Dirac [27| . and since then many papers have appeared 
on this subject j28l - l32| . Nevertheless, there are only few 
classical electrodynamical studies on the standing wave 
ponderomotive force ^33,-^]. Most publications on the 
ponderomotive force have concentrated on propagating 
EM waves, establishing the validity of Eq. ^ in that 
context [l|, 0, [l3| • Eq. ((TJ is also applied to the standing 
wave configuration, implicitly assuming that it remains 
valid in that case as well. In 2005, however, Kaplan 
and Pokrovsky [s^, [11] calculated the time-averaged 
equation of motion of an electron in a standing wave 
for a number of field polarizations, and their results 
showed that the ponderomotive force depends on the 
polarization. Most notably, the ponderomotive force can 
even change its direction toward high field regions for 
certain situations. Clearly, these results are in confiict 
with the polarization-independent Eq. ^ that is 
commonly used. Kaplan and Pokrovsky did not provide 
an alternative expression for the ponderomotive force, 
however. 

We would now first like to show, on the basis of 
simple arguments, that it can be understood that the 



In this paper, the ponderomotive force is derived for 
a relativistic particle entering a nonrelativistic standing 
wave with a general three-dimensional field distribution. 
It is shown that indeed Eq. ([T|), or the relativistic equiv- 
alent, does not hold for this situation. This may have im- 
portant implications for experiments and proposals based 
on the standing wave ponderomotive force. The main re- 
sult of this paper, Eq. (HJ), shows that the modified 
ponderomotive force is still of simple gradient form, but 
contains additional polarization-dependent terms. We 
identify the polarization-dependent contribution due to 
the magnetic quiver. An additional relativistic effect ef- 
fectively doubles this magnetic contribution. We thus 
generalize the results of Kaplan and Pokrovsky, which 
follow naturally from our ponderomotive force expres- 
sion. This paper is structured as follows. In Section 
im the polarization-dependent ponderomotive gradient 
force in a standing wave, Eq (j21[) . is derived using a 
perturbation expansion method. The origin of the addi- 
tional polarization-dependent terms is discussed. Next, 
Eq. (PT|) is applied in Section IIIII to calculate the aver- 
aged equation of motion of an electron in specific stand- 
ing wave geometries, reproducing the results obtained by 
Kaplan and Pokrovsky. In Section IIV[ we validate our 
ponderomotive force expression by testing it against nu- 
merical simulations of electron trajectories in a realistic, 
non-idealized field configuration. 
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II. THE POLARIZATION-DEPENDENT 
PONDEROMOTIVE FORCE 

A. Assumptions 

First consider the idealized case of two plane EM waves 
of equal frequency to that counterpropagate along the 
2-axis and add to form a standing wave. In the Coulomb 
gauge, the vector potential A of this ideal standing wave 
is then purely harmonic in time t and position z, i.e. 



A = and [p'^ jdz^ + /c^) A — 0, where 



k = Lo / c with c the speed of light. Furthermore, the 
vector potential satisfies — and V ^_Ai = 0, where 
i = x,y. 

In practical applications, however, an EM standing 
wave differs from this idealized situation in two ways. 
First, the standing wave has a finite transverse extent, 
leading to a small but nonzero transverse gradient V ±Ai 
and a small longitudinal component A^ . These two quan- 
tities are related by the gauge condition V • A = 0, and 
scale analysis of the latter shows that, symbolically, 



A^ 

Ai_ 



d/dz 



(2) 



Here, e ^ 1 is a small parameter measuring the magni- 
tude of the field inhomogeneity, and will be used as the 
expansion parameter in the derivation that follows. For 
example, in Gaussian laser beams focused to a waist of 
size wq, this parameter is e ~ (fcwo)^^- 

Secondly, the standing wave has a finite lifetime, so 
that the vector potential is only quasi-monochromatic: 



OiSA); 



d^A 

d{kzy 



0{SA), (3) 



where ^ 1 is another small parameter measuring the 
monochromaticity. For a standing wave produced by 
counterpropagating laser pulses of temporal length a, 
for example, this parameter is 5 ~ (wtr)"^. 

In addition, in this paper the EM field is assumed to 
be of nonrelativistic intensity, which means that 



eA 



mc 



a < 1. 



(4) 



For reasons of clarity, for the moment it is assumed that 
a ~ e. However, the derivation below can be generalized 
straightforwardly to other field strengths such as a ~ 
or a ~ e^^^, leading to the same result. Appendix IB] gives 
a short description of the generalized derivation. 



B. Perturbation expansions 

Before considering the dynamics of a charged particle 
in a standing wave in detail, let us first determine what 



time scales are involved. First, there is the time scale of 
the quiver motion, which is the optical time scale 
Second, there is the time scale on which the motion of 
the guiding center changes. Substitution of an ideal 
standing wave A — AqBx cos kz sin ujt in Eq. ([T]) and 
integrating yields oscillatory motion in the z-direction, 
with a typical frequency = cAquj / {\/2mc) . Thus the 
guiding center motion in the z-direction changes on a 
second, longer time scale (auj)^^ ^ (euj)^^. Finally, 
in a realistic standing wave, the nonzero transverse 
field gradient causes transverse ponderomotive forces, 
which in view of Eqs. ([T]) and ^ are weaker than 
the longitudinal ponderomotive forces by a factor e. 
Therefore the transverse guiding center motion changes 
on a third, still longer time scale (e^w)~^. 

Having established the three time scales of the prob- 
lem, next consider the equations of motion of a charged 
particle in the standing wave ^39] : 



(5) 
(6) 



d{ijt) 



(-- 




1 /AVeA\ 


P 


\mc 


mc ) 


7 \ 27r mc J 


mc 




d{kx) 








d{ujt) 


7 mc' 





in which 7 = -^/l + p'^/ (mc)^ is the Lorentz factor, 
and the dyadic notation VA has been used [40]. Eqs. 
(ini)-® have been made dimensionless by dividing the 
usual equations by mcuj and c respectively. Below, these 
equations are solved by expressing the various quantities 
in perturbation expansions in terms of e. Subsequently, 
terms of like order in e will be collected and equated 
(4ll |. We use the symbol 'Ord' to denote 'on the order 
of, i.e. / = Ord{e^) means < limej^oZ/e* < 00, 
in distinction with / = 0(e*) which is equivalent to 
< linijj^o //f * < 00 ■ Superscripts in parentheses denote 
the order of the terms. 



First, the momentum is expanded as 



P 



P 



(i) 



1=0 

Ord{e'). 



(7) 



Note that would be the momentum in the absence 
of the standing wave field, that is, the initial momentum, 
since p — > p'"-* as a ^ e | 0. Next, each order p^*' 
is decomposed into a slowly varying guiding center 
part p*^*^ — (p^*'') and a rapidly varying quiver part 



pyi — py/ — [p'-^'). Here, (•) denotes time-averaging 
on the time scale u!~^. Upon substitution of this 
decomposition in the left-hand side of Eq. ^ , each term 
is differentiated with respect to ut. This preserves the 
order of magnitude of the fast quantities p^*) , since these 
vary on the time scale w"^, i.e. dp'^^'' / d{u>t) — Ord{e^). 
However, from the discussion above, the slow quan- 
tities pi*' and p^-* vary over longer time scales, so 
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that differentiation increases their order according to shown below. 



dpr/d{ujt) = Ord(e*+i) and dp^l> /d{ujt) = Ord{€'+^) 
respectively. More formally, these order relations may 
be established using the multiple scale technique [4l|, 
considering pz a function of the two variables ut and 
ecut and considering p± a function of the two variables 
ujt and e^Lut. Multiple scale analysis has been applied 
in a relativistic derivation of the ponderomotive force in 
propagating EM radiation (T^ . 



Next, the expansion ([7]) is substituted in the recipro- 
cal Lorentz factor I/7. Extracting the zeroth order part 

(0) ^ 



7 



1 



1 

7 



^p(o)j /(toc) 2, this gives 
1 



7(0)^1 + (mc7(o)) (2p(")-p(i) 



7 



(0) 

7) + 



p(o) ._p(l) 
(mc)2 (7(0) 
1 



(1) 

-] +0{e'). 



(8) 



The quantity (l/7)'-^-' is the first-order time-dependent 
variation of the reciprocal Lorentz factor with respect to 
the constant value (1/7)^°^. As will be shown below, the 
associated variation in the relativistic mass 7m of the 
particle leads to an additional contribution in the final 
ponderomotive force expression. 

Finally, below it will be required to take the time av- 
erage of expressions involving powers of A or its deriva- 
tives. These time averages need to be taken along the tra- 
jectory of the particle; that is, in the average {A {x{t), t)), 
the vector potential is to be evaluated at x — x{t). To 
bring out this position dependence explicitly, the position 
is also expanded in a perturbation expansion in terms of 
e: 



i=0 



(9) 



= Ord (p(*') = Ord{^) 



dt 



Again, each order a;(*) is decomposed into a slowly vary- 
ing guiding center part S^'-* — {x'^'^^) and a rapidly vary- 
ing quiver part a;'^') — x^^^ — {x^^'). Then the vector 
potential can be expanded in a Taylor series around 

X 



dA 



A{x{t),t) = A + I^'\t)^ + Oie') 

oz 



(10) 



Here and in the remainder of the paper, an overbar on the 
vector potential denotes evaluation at the guiding center 
position, i.e. A = A{x,t). In writing the series in Eq. 
(ITU)) , it has been anticipated that ir^"^ = 0, as will be 



C. Order-by-order solution of equation of motion 

Substitution of the expansions ([71)- (|10l) in the equa- 
tions of motion (HI)-®, and collecting terms of equal or- 
der in e, results in two equations at each order of e. These 
order equations are listed in Appendix[Xl Order-by-order 
solution, balancing in each equation the averaged parts 
and the oscillating parts separately, yields the zeroth or- 
der quantities 



P 



(0) 



•(0) _ 



P 

dx^ 
dt 



0: 



Po 



7717 





(0)' 



(11) 
(12) 

(13) 

(14) 



in which pq is the initial momentum. As expected, at 
zeroth order (that is, in the limit a ^ e ^ where both 
the field strength and field inhomogeneity are zero) the 
motion is equal to what it would be if the EM field were 
absent. For the first order quantities, it is found that 



(15) 

■Po±; (16) 



p'l" 


= -eA± 




dpi 


e 


dAi_ 


dt 


7717(0) 


dz 




1 




dt 


7777(0) 








dt 







-POz 



eAj_-poj 



PazPz 



(mcf 



(0)1 



;(i7) 

(18) 



Eq. (jl5p expresses the well-known result that, in an os- 
cillating EM field, at lowest order the quiver momentum 
balances the vector potential, such that the canonical 
momentum p + eA is conserved. Eq. (fT5)) will be used 
in Section IIIII for the description of the guiding center 
motion. Eqs. (|16p - (|17p describe the quiver motion in 
the direction normal to the plane of polarization of the 
standing wave, which is the direction of strong field gra- 
dient. This is the magnetic quiver motion described in 
the Introduction. Time differentiation of Eq. (flT)) and 
substitution of Eq. ([TO]) yields 



dt^ 



(to7("))2 



(l-/3o.)^ + 
oz 



dA±_ , /3o^ dA^ 



dt 



■Po±, 



(19) 



where /3o = Po/i'mc^^^^) is the initial velocity divided by 
c. We now restrict to the situation that /3oz is sufficiently 
small so that the second term in parentheses in Eq. (fTU)) 
is negligible, which is the case if /3oz ^ 1- Then, in 
addition, using Eq. ([3]) the full time derivative may be 
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written (P/dt^ = {d/dt + v-Vf = -uj^ + 0{e,S, (Sq^), so 
that double time integration yields 



j(i) = _ 



Po±[l + 0{6,Poz)]- (20) 



The bracketed factor expresses the error introduced 
by the integration. Eq. ^IU\\ clearly shows that the 
amplitude of the magnetic quiver motion is polarization- 
dependent: if Po± is parallel to A±, this amplitude may 
be substantial, while for po± perpendicular to A_l, it 
vanishes at first order. This reinstates the argument 
made in the Introduction: if po_L |1 A±, the momentum 
is largely perpendicular to the magnetic field, resulting 
in a substantial magnetic force and quiver amplitude. 
Conversely, if poJ. -L the momentum is largely 

parallel to the magnetic field, with vanishing magnetic 
force and quiver amplitude. 

As a final step, we substitute Eqs. (fTT|) - (IT5|) and ((20|) 
in the right-hand sides of the remaining order equations 
(IA7|) - (IA8|) . and take the time average of these equations. 
Then the left-hand sides reduce to the rate of change of 
the first order guiding center momentum, dp^^^ /dt. The 
right-hand sides reduce to a single gradient: 



dt 



(21) 



•V(a1-(/3, 



•0_L 



■A, 




This is the polarization-dependent ponderomotive force 
in a nonrelativistic standing wave for a particle with 
Poz ^ 1; it is the main result of this paper. The approxi- 
mate sign expresses a relative error of the order of /3oz + 5. 
In the limit (3o± 0, Eq. (PT|) reduces to the well-known 
polarization- independent ponderomotive force, Eq. ([ij, 
with the relativistic factor 7^°) « \/l -f Jp / (mc)^ + a^/2 
included in the denominator. For /3o± 7^ 0, however, the 



two polarization-dependent terms of Eq. 
significant. 



211) become 



The term in Eq. ([2T|) proportional to 

(9 (/3oJ_ • Aj_) /dkz) originates from including the 
magnetic quiver motion, Eq. (j20p . in the Taylor ex- 
pansion of the vector potential, Eq. (ITU]). It accounts 
for the fact that the z-position oscillates in phase with 
the temporal oscillation of the field. Therefore, when a 
field gradient in the z-direction is present, the particle 
systematically samples higher fields at selected phases 
of the electric quiver motion in the direction of A^, 
and lower fields at other phases. The induced nonzero 
average force is negligible in most applications of the 
ponderomotive force. But in a standing wave the 
derivative in d'^/d{kzY is of order unity, so that this 
term is comparable to the other terms in Eq. (j2ip . 



The origin of the term in Eq. ([2T|) proportional to 
(/3ox • A±) can be traced back to including the first or- 
der term (1/7)*'^^ in the expansion of the Lorentz factor, 
Eq. ([S]), rather than approximating 7 « 7'^'^) throughout. 
This takes into account that the relativistic mass 7m of 
the particle oscillates around the value ^'''^^m in phase 
with the oscillation of the field. Therefore, during its 
primary quiver motion in the direction of Aj^, the par- 
ticle is systematically heavier at selected phases of the 
quiver motion, and lighter at other phases, leading to a 
nonzero effective force when a field gradient is present. 



III. WIGGLING MOTION IN A STANDING 
WAVE 



The last term of the z-component of Eq. (|2T|) may be 
rewritten by performing the z-differentiations and using 
Eq. ([3]), after which Eq. (PT|) becomes 



dp. 



(1) 



dt 



2m7(0) dz 



(22) 



We now evaluate Eq. (l22t for the linearly and circularly 
polarized standing waves considered by Pokrovsky and 
Kaplan. 



A. Linear polarization 

Let the standing wave be produced by two counter- 
propagating plane waves of equal amplitude and fre- 
quency that are coUinearly polarized in the direction 
Cp — exC0s4' + eySvucj), as is indicated in Fig. [T] by 
the dashed arrow. Then the vector potential is A = 
AqBp cos kz sin ojt. Suppose that a charged particle en- 
ters the standing wave parallel to the x-axis with initial 
velocity /Jg. Differentiating Eq. (|18l) . and substituting 
Eq. ([m, gives the equation of motion for the guiding 
center of the particle in the z-direction: 



diojtf 



2 1 - 2/3^ cos2 (j) 
'° 4(7("))' 



sin(2fcz) « 0, (23) 



in which flo = eAo/(mc). Eq. ((23)) shows that the guiding 
center makes pendulum-like oscillations in the z-direction 
(it wiggles in the ponderomotive potential), with equi- 
librium points at fcz = n7r/2 and a small-amplitude 
frequency = aoUj^\l - 2/3^ cos2 0|/(\/27(o)). 

First, consider polarization parallel to the initial veloc- 
ity, i.e. cos(/) = 1. Then Eq. ^ reduces to Eq. (40) of 
Ref. [11] after rewriting /3o = po/(wc7^°^). As was noted 
in Ref. [11], the most striking feature of this configu- 
ration is that the guiding center oscillations vanish for 
Pa — 1/ v^; in terms of the guiding center motion, the 
standing wave is invisible to the particle for this value 
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of initial velocity. When /3o is increased above l/-\/2, the 
stable and unstable equilibrium points of Eq. ()23|) reverse 
their positions, i.e. the ponderomotive force changes di- 
rection towards high field regions. This is the relativistic 
reversal described in Ref. [SSj. 

For polarization perpendicular to the initial velocity of 
the particle, cos(/) — 0. Then Eq. ([23l) becomes identi- 
cal to Eq. (52) of Ref. [s^. For this polarization, the 
magnitude of the ponderomotive force and the wiggling 
frequency are independent of the initial velocity, and the 
relativistic reversal effect is absent. 



B. Circular polarization 

If the standing wave is produced by two counter- 
propagating plane waves that are circularly polarized 
with opposite helicities, the vector potential is equal to 
A — Ao cos kz (fix coswt -I- Cj, sincji). Thus the standing 
wave has equally spaced nodes and antinodes along the 
z-axis, while locally the field direction rotates around the 
z-axis with time. Again combining Eqs. (|T8| and ((22)) 
and substituting the vector potential now gives 



d(ujt) 



v(o)l 



■ sin(2fcz) « 0, 



(24) 



which equals Eq. (56) of Ref. [38)] taking into account 
that the field amplitude used there is AqI\/2. Thus, in 
this configuration the equilibrium points are the same as 
in the linearly polarized case described by Eq. for 
cos </> = 0, although the magnitude of the ponderomotive 
force is a factor (7''"-') weaker. 

For counterpropagating circularly polarized 
waves of equal helicity, the vector potential reads 
sin (e^; sinfcz -I- By cos fcz). In this case, the stand- 
ing wave is polarized in a helix along the z-axis. Eqs. 
(UHl) and ^ yield 



d(ujtf 2 (7(0)) 



- sin(2fcz) w 0, 



(25) 



which is the same as Eq. (62) of Ref. [3^ after rewriting 
/3o ~ Po/("^c7*-'^-') and taking into account the different 
definition of the field amplitude. As was noted in Ref. 
[ssl ]. the field intensity is homogeneous along the z-axis, 
so that Eq. ([T]) predicts zero ponderomotive force. 
However, the modified ponderomotive force expression, 
Eq. (pT|) . shows that this force is nonzero, so that Eq. 
(|25p still yields wiggling motion of the guiding center. 



In summary, all the equations of motion (P ^ - (P5)) can 
be derived from a single ponderomotive force expression, 
Eq. ([21]) . We have therefore generalized the results of 
Kaplan and Pokrovsky, who started from the Lorentz 
force equation for each individual case. Moreover, Eq. 



([21]) also gives the transverse component of ponderomo- 
tive force, which was not considered in Ref. [38|. 



IV. COMPARISON WITH SIMULATIONS 



In this section, Eq. ([2T|) is tested against numerical 
simulations of electron trajectories in a realistic, non- 
idealized standing wave field. We have used the GPT 
code which uses an embedded fifth order Runge-Kutta 
method with adaptive stepsize control 42? ]. For 
comparison with the simulations, Eq. (j21l) is needed 
in terms of the electric field rather than the vector 

potential. From Eq. ([3]), the potential is approximately 

2 2 

harmonic in time, so that (A^) « lli^'^{E^). Therefore 

Aj_ may be effectively replaced by E^/lo in Eq. (|2l]) . 

We consider again the configuration shown in Fig. 
[1] this time with two identical Gaussian laser beams 
in the fundamental mode that counterpropagate along 
the z-axis, have a central wavelength A = 800 nm 
and a peak intensity /q = 2.0 • 10^^ W/cm^, and are 
focused in a common waist of size wq = 12.5 /im at 
z = 0. For these parameters, e ~ {kwo)^^ — 0.01 and 
a ~ 0.01 so that the theory of Section [TTl is valid. The 
beams are assumed to be pulsed with a Gaussian pulse 
shape of length ct = 85 fs, and timed such that the 
centers of the pulses coincide at z = at time i = 0. 
This configuration yields a standing wave around the 
origin with length ~ ccr and width ^ iuq, and with a 
lifetime ~ a around t = 0. Furthermore, we assume 
that the pulses arc coUinearly polarized in the direction 
Bp — Ex cos (p + By sin 0, similar to the case considered 
in Section IIII Al In each of the following cases, the 
trajectory is calculated of an electron entering with 
initial velocity /3o = [ioex for t <C — cr, and with the 
initial position chosen such that, at i = 0, the position 
of the electron would be {x,y,z) = (0, Wo/2, A/8) in 
absence of the laser fields. Then the electron meets 
the standing wave close to the origin, has interaction 
with it for some time ~ a, and leaves the interaction 
region in a defiected direction. The initial position has 
been chosen such that the electron samples the highest 
available field gradients both in the z-direction and in 
the perpendicular direction, in order to maximize the 
ponderomotive effects. 

In Figs. [2]and[31 the trajectory of the electron is shown 
as viewed from the negative y-axis, for several values of 
/?o and for two different polarization directions. The tra- 
jectories are calculated using three different methods: 

• (GPT) numerical integration of the exact equations 
of motion using the GPT code (solid lines); 

• (OLD) calculation according to the classical pon- 
deromotive force, Eq. ([T]) (dashed lines); 

• (NEW) calculation according to our polarization- 
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dependent ponderomotive force, Eq. ([21 
replaced by Ej^juj (dash-dotted lines). 



with A I 



Standard paraxial field expressions have been used [l8l |: 
these are listed in Appendix [C] for reference. 

Fig. [2] shows the configuration in which the polar- 
ization is perpendicular to the initial velocity. In this 
case the polarization-dependent terms of Eq. (|21|) van- 
ish, and Eq. (PT|) reduces to the classical expression, Eq. 
([1]). Indeed, for all three initial velocities the three de- 
scriptions yield identical trajectories to within the width 
of the lines, showing that the classical expression gives 
an excellent description of the averaged motion of the 
electron. The usual behavior can be seen in which the 
electron is defiected towards low-intensity regions. 

In Fig. [31 however, the situation is very different. Here, 
the polarization is parallel to the initial velocity, so that 
the polarization-dependent terms of Eq. ([2T|) become 
important. Because of this, in Fig. [3l^a) the magnitude 
of the ponderomotive force is smaller than in the corre- 
sponding case of Fig. [2] (/3o = 0.5), the ponderomotive 
force vanishes in Fig. [Sjb), and it even changes direction 
towards the high- intensity region in Fig. [Hfc). In all of 
these cases, the resulting trajectories are excellently pre- 
dicted by Eq. (|2ip . Fig. [3] de monstrates the relativistic 
reversal described in Ref. [3g| and Section IIII Al Mean- 
while, Eq. ([!]) is insensitive to the polarization direction, 
so that the trajectories (dashed lines) are incorrectly pre- 
dicted to be identical to the corresponding cases in Fig. 

m 

The insets of Fig. [31 are close-ups of the solid lines, 
showing that these actually consist of the smooth, time- 
averaged trajectory predicted by Eq. ([2T1) (dash-dotted 
line), and the GPT trajectory. The latter contains 
also the quiver motion, which has components in the 
x-direction (Eq. (ITS]) ') and the z-direction (Eq. ([^ ) 
and is therefore visible in the plane of drawing of Fig. 
[3l The frequency of this quiver motion is uj while the 
electron moves forward with a velocity /3x ~ /^o^i so 
that the spatial period of the quiver motion is l/(/3oa;A). 
This has been made visible by scaling the x-axes of the 
insets. In Fig. [21 the quiver motion would not be visible 
in a close-up because it is perpendicular to the plane 
of drawing, since ■ po_L = in Eq. ([20l) and the 
z-component of the quiver motion vanishes. 

In order to show some ponderomotive effects in the 
perpendicular direction as well, in Fig. [4l a sideview is 
given of the trajectories of Fig. [SJa), as seen from the 
positive z-axis. Note from the vertical scale that in this 
direction the defiection of the electron is very small due 
to the very small field gradient. Nevertheless, again it is 
clear from the figure that Eq. ([2T|) accurately predicts 
the electron trajectory, contrary to the classical expres- 
sion. Thus Eq. ([^ gives a precise description of the 
three-dimensional electron trajectory in a realistic, non- 
idealized field configuration. We have also repeated the 
simulations for combinations of e and a other than a ~ e. 
Whenever both e and a are less than about 0.1, consis- 



tent with the assumptions a, e ^ 1 made in Section[TTl we 
find the same level of agreement with the GPT results. 



V. CONCLUSION 

The classical polarization-independent ponderomo- 
tive force expression is commonly used to describe 
the time-averaged motion of a charged particle in an 
inhomogeneous oscillating EM field. It is generally 
assumed that this is an accurate description, at least for 
nonrelativistic field intensities. However, we have shown 
that this is not always true. If the field configuration 
possesses a direction in which the field changes on the 
scale of the wavelength, i.e. in a standing wave, and 
in addition the charged particle is relativistic, the pon- 
deromotive force is modified. In particular, it becomes 
dependent on the polarization of the field. Because 
of this, the ponderomotive force may even vanish, or 
change its direction toward high field regions, as was 
found earlier by Kaplan and Pokrovsky [37]. We have 
derived the modified ponderomotive force expression for 
these configurations, which is of gradient form like the 
classical expression. Comparison with simulations in 
the case of a realistic, non-idealized, three-dimensional 
field configuration confirmed the general validity of the 
analytical results. 

The modifications of the ponderomotive force derived 
in this paper may have important implications for ap- 
plications that involve the ponderomotive interaction of 
relativistic charged particles and standing EM waves. 
For example, in the electron bunch length measurement 
based on ponderomotive scattering of the electrons by a 
standing wave [19], the polarization of the wave is es- 
sential for an optimal design of the experimental setup. 
In the proposed X-ray free electron laser relying on the 
wiggling motion of electrons induced by the ponderomo- 
tive force in a standing wave |2l| , the frequency of wig- 
gling and hence that of the stimulated radiation directly 
depend on the polarization. Experimental tests involv- 
ing the controlled scattering of electrons by a standing 
wave have confirmed the classical ponderomotive force 
expression |28i] and Kapitza-Dirac diffraction [29] using 
nonrelativistic electrons. It would be very interesting 
to extend these experiments to relativistic electrons to 
test the polarization-dependent ponderomotive force, Eq. 
([^T|). in the classical limit, and to study the polarization- 
dependence of Kapitza-Dirac diffraction. 



Appendix A: Order equations 

Substitution of the expansions ©-([TU]) in the 
equations of motion ([S])-®, and collecting terms 
of equal order in e, results in the following order 
equations. The components of the vector poten- 
tial and the spatial derivatives have been treated 
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as ej4j^/(mc) = Ord{e), eAz/{mc) = Ord{e^), clarity, the equations are displayed in dimensional form. 
AVx = Ord{e) and Xd/dz — Ord{l) consistent 
with Eq. ([2]) and the assumption a ^ e. For reasons of 



Ord{l) 



Ord{e) 



dp(°^ 

dt 


= 0; 






d{x + x)^°^ 


{P + P)f 






dt 








dp^}'^ dA±_ 
-^-^ + e 


= 0; 






dt dt 








e dA±_ 


•(P- 




dt 




d{x + xf^^ 






^ p)^"^ ■ {p + pf^ {p+pf^ 


dt 









dA, 



Or die") 



dt 



(0) 



(V^Ax)-(p + p)f ; 



(0) 



{P + P)j 



{p + P)'^ + ^-^{p + P)T- 

■{p+P)T ■ 

{VaA^) ■ {p + p)^l^ + (V^A,) (p + p) 

(0) 



dz^ ' ^'^ ' dz 



(mc7(°))^ 



9z 



(0) 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 



dz 



(A7) 



(0) 



(v^Ax)-(p+p)r 



(A8) 



Appendix B: General field strength 

The ponderomotive force, Eq. (f2T|) . has been de- 
rived under the assumption that a ^ e. Since this 
expression is the result of balancing terms of equal 
order in the equations of motion, one might expect 
that it would be affected by changing the order of 
magnitude of the vector potential to, for instance, 
a ~ e^. This is not the case, however. We only give 
a sketch of the generalized derivation for arbitrary a 1. 

Repeating first the order expansion method for the 
case a ^ e"", n > 1, it is not difficult to find the 
lowest order slowly varying term and the lowest order 
rapidly varying term of the expansions of p and I/7. 



Also, the first two terms of the expansion of A follow 
straightforwardly. The right-hand side of the equation 
of motion ([5]) is then formed by factoring out the 
product of these three expansions. Taking the time 
average of the result, it is found that the lowest order 
terms that are nonzero on average are precisely those 
that form the ponderomotive force given by Eq. pip . 
In terms of the corresponding set of order equations 
analogues to those in Appendix \^ the first 2n — 1 
orders of the momentum equations yield zero right-hand 
sides upon averaging, while the 2n-th and (2n + l)-th 
orders evaluate to respectively the z-component and 
perpendicular component of Eq. (j21l) . 

Also the opposite situation in which 1 ^ a ^ e^/". 
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n > 1, is possible. Since in this case factors of A in 
the equations of motion lead to terms of fractional or- 
der in e, it is appropriate to expand all quantities in 
power series in terms of e^/" rather than e. This is ef- 
fected by using the same power series expansions as be- 
fore, with the understanding that p^*^ = Or(i(e'/") rather 
than = Ord(e'), for example. Except for this modifi- 
cation, the derivation of the ponderomotive force is anal- 
ogous to that for the case a ~ e" considered above, and 
again Eq. (PT|) is found. Thus Eq. (PT|) is valid for arbi- 
trary a ^ 1. 



Appendix C: Field expressions used in numerical 
calculations 



The solid lines in Figs. [2M] have been calculated by 
numerical integration of the equations of motion 



in which Eq = ^ 2Iq / [cqc) is the peak electric field ampli- 
tude, w = wq^J1 + I z\ is the beam waist, zr = kuy^jl 
is the Rayleigh length, and the Gouy phases are 



X± = bJt±[k 



arctan h 

z 

2 arctan — 



z X + 
zr 

z + y"^ 
H — 

Zr 



The dashed lines in Figs. have been calculated ac- 

2 

cording to Eq. ([!]) with I{x, t) = eac{E^) using Eq. (jClj) 
for Ex- The dash-dotted lines have been calculated ac- 
cording to Eq. (|2ip with A± replaced by Ex/to. For 
polarization in the y-direction (0 — 7r/2), replace x ^ y 
and y ~x in the expressions above. 



^ = e ( E + -p X B 
dt \ 1 

dx p 

dt mj ' 



using for E and B the following paraxial Gaussian beam 
fields For polarization in the x-direction (0 = 0), 



E = £;+ exp - 



B = B+cxp[- 



{z + ctf 
{z + ct)^ 



E^ exp — 



B cxp — 



jz - cty 

(CI 

{z - ctf 



E± = Eq — Bx cos V'± ± smx± ; 

w \ ZrW 

n Eow„ f ywo 

B± = [^By cos ^p± e^smxi 

C W \ ZrW 
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Polarization out of plane (p = nil) 




^ ^ ^ I 1=0 

FIG. 2: Trajectory of an electron incident on a standing wave 
that is polarized in the j/-direction (out of plane), for three 
different initial velocities /?o- The color map shows the field 
intensity of the standing wave in the {z, a;)-plane at time t = Q. 
For each initial velocity, the plot actually includes three tra- 
jectories, calculated with the methods GPT, OLD, and NEW 
respectively; in each case all three lines overlap to within the 
linewidth. 




FIG. 3: Trajectory of an electron incident on a standing wave 
that is polarized in the a;-direction (left-right), for an initial 
velocity /3o equal to (a) 0.5; (b) 0.707 « 1/^2; (c) 0.9. The 
color map shows the field intensity of the standing wave in 
the {z, a;)-plane at time t — 0. In each plot, the trajectory 
has been calculated with the methods GPT (solid line) , OLD 
(dashed line) and NEW (dash-dotted line). 
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FIG. 4: Sideview from the positive 2:-axis of the trajectories 
shown in Fig. [3ja). The trajectories has been calculated with 
the methods GET (solid line), OLD (dashed line) and NEW 
(overlapping with GPT). 



